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ABSTRACT 



We investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. 
To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and 
different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is 
a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and 
mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased 
by a factor larger than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are 
also rather short-lived, with a typical duration of the activity of a few 10* yr. These conclusions are found to be robust, independent 
from the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced 
starbursts are neither stronger neither longer than their local counterparts. In turn, the formation of massive gas clumps, results of 
local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more 
important role in determining the duration and intensity of star formation episodes. 



1. Introduction 

The role played by ga laxy interactions i n affe cting star for- 
mation was realized by iLarson & Tinsleyl (Il978[). who sh owed 
that disturbed galaxies in the Arp Catalogue (lArp|[r966l) have 
a larger dispersion in their colors and a bluer envelope in the 
, (U-B, B-V) plane than normal systems taken from the Hubble 
• Atlas dSandagd Il96lh . Using evolutionary synthesis models, 
' they suggested that the features found in U BV colors of inter- 
acting galaxies were caused by bursts of star formation lasting a 
few 10^ - 10** ye ars. The large am ount of observational works 
that followed (see iKennicutt et al] (11996) for a complete review) 
, showed that in the Local Universe, the response of galaxies to 
mutual interactions and mergers is quite varied. 

Many starbursts in the Local Universe take place in the cen- 
tral regions of interacting/merging gal axies, as it is the cas e 
for instance for NGC 7714 studied by IWeedman et al.1 (1 1981 1 ) . 
and the protoype starburst galaxy M 82 (Ide GriisI l200Tallbl) . 
Another well-studied example is the NGC 4038/4039 system 
(the Antennae): this early stage merger presents an extended star 
formation, the most intens e star forming r egions being located 
between the two galaxies (IWang et al.ll2004 ). Actually, the vast 
majority of UltraLuminous Infrared Galaxies (ULIRGs) at low 
redshift, i.e. the strongest starbursts in the Local Universe, are 
found in interac ting and m erging galaxies (Sanders & Mirabel 
(119961) . see also iDuc et al . I (119971)) This is however not recip- 
rocal. Indeed, iBergvall et al.l (l20'04b among others have shown 
that, in a magnitude-limited sample of 59 interacting and merg- 
ing galaxies, only a weak enhancement of star formation (a fac- 
tor of 2-3 in the galaxy centers) is found when compared to a 
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reference sample of non-interacting galaxies, so that the contri- 
bution of interactions and mergers to the global star formation 
activity at low redshift could on average be much less efficient 
than suggested by the strongest examples of starbursts. 

At high redshift, the role of mergers in the star formation his- 
tory of galaxies i s debated, too. In a pioneering study of distant 
infrared galaxies. ^ Elbaz & Cesarskvl (|2003[) showed that the ma- 
jority of present-day stars were formed in dusty starbursts, and 
suggested that t he late r were triggered by galaxy interactions. 
IConselice et al.l (l2003l) also suggested that about two thirds of 
submillimeter galaxies at z > 1 are undergoing a major merger . 
In a study of the Spitzer First Look Survey, Bridge et al. (2007|) 
argued that close pairs are major contributors to the star forma- 
tion density at z > 0.7, about half of the star formation rate den- 
sity at z ~ 1 being attributed to major mergers and interactions. 
However, an important part of the infrare d luminosity of g alax- 
ies could be caus ed by AGN heating (e.g. lDaddi et alj|2007allbl) . 
iBell et alJ (l2005h find that less than one third of actively star 
forming galaxies at z ~ 0.8 are actually interacting or merging, 
the majority having undisturbed disk morphologies. Similarly, 
selecting interacting galaxies in the GEMS survey, Jogee et al] 
(l2007l 12008 in preparation!) find that over the redshift interval 
z~0.24 to 0.80 (corresponding to lookback times 3 to 7 Gyr), 
the average SFR of strongly distorted interacting/merging mas- 
sive galaxies is only modestly enhanced with respect to nor- 
mal undisturb ed galaxies. At even higher redshift (z ~ 2) in 
GOODS, Da ddi et a n (l2007al) find that the star formation activ- 
ity of ULIRGS is long lived (at least half a Gyr) which might be 
longer than expected for merger-induced starbursts. 

In the last decade, a lot of progress has been made in our 
theoretical understanding of the role played by galaxy interac- 
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tions in driving star fo rmation, and a large variety of results 
has been obtained (see IStruckl (|2006|) for a review). Barnes & 
Hernquist (1991) have shown that tidal interactions between 
galaxies can drive gas inflows toward s the central regions, which 
can in crease the star formation rate. Mihos & Hernquisl (Il994ai 
1 19961) directly studied the star formation activity in mergers of 
equal mass disk galaxies using self-consistent N-body simula- 
tions (jMihos & Hernquist 1994b) including stars, dark matter, 
and gas dynamics. They also pointed out that the timing and 
strength of an interaction-driven starburst depend on the mor- 
phology of the interacting systems, bulgeless disk galaxies being 
more prone to suffer enhanced star formation in the first phases 
of the interaction, rather than in the merging phase. Important 
episodes of sta r formation can also arise in the case of minor 
merger events dMihos & Hernquist 1 1994cl) . but major mergers 
of galaxies of comparable masses are t he most efficien t situa- 
tions to trigger strong starbursts. Indeed. ICox et al. I ('2008) have 
shown that the star formation activity of merging galaxies de- 
creases rapidly with increasing mass ratios. More sophisticated 
models including supernovae feedback in regulating star forma- 
tion have also been explore d (Springel 2000; Cox et al. 2006, 
12008 ). In particular. [Cox et al.1 (12006) pointed out that the large 
amount of freedom in selecting the feedback parameters could 
play a significant role in determining the maximum star forma- 
tion rate during a galaxy merger Indeed, supernovae feedback 
regulates the star formation and, on average, reduces the inten- 
sity of merger-induced starbursts. The effect of feedback is how- 
ever generally modest within the assumptions considered as the 
most realistic ones. 

Nevertheless, that some cases of major interactions or merg- 
ers can trigger strong star bursts, as in the examples shown by 
iMihos & Hernquisa(ll996l) . does not imply that the enhancement 
of the star formation activity is systematically high. A sample 
of about 50 Nbo dy-SPH simu l ations of galaxy interactions was 
performed by IKapferer et al.l ( l2005l) . and these authors found 
that the integrated star formation rate during an interaction is 
moderately increased, up to a factor of 5 but on average a fac- 
tor of 2 with r espect to that of isolated galaxies. More recently, 
iDi Matteo eT al. (2007) (hereafter DM07) presented more than 
two hundreds simulations of galaxy interactions and mergers, 
restricted to coplanar cases, pointing out the difficulty to drive 
intense starbursts. Only 17% of the mergers in their sample have 
strong bursts with an SFR ten times higher than in isolated galax- 
ies, and half of the sample shows an enhancement of the SFR by 
a factor no bigger than 4 at the peak of the starburst. 

In this paper, we extend the analysis started in DM07 on the 
relation between major galaxy interactions and star formation. 
On the one side, we enlarge the sample of Tree-SPH simulations 
studied in that paper, modeling orbits with various inclinations. 
On the other side, we test whether the main conclusions depend 
on the model used. To this aim, we compare the result of this first 
sample to a second (smaller) sample of simulations performed 
with a different numerical code - a particle-mesh sticky-particle 
(hereafter PM-SP) - using somewhat different initial conditions 
and testing different numerical recipes for star formation. This 
large and heterogeneous data set shoud help in finding robust 
results about the starbursts-galaxy interaction connection. It is 
also, to our knowledge, the first numerical work in this field 
where simulations realized with different codes and star forma- 
tion recipes are directly compared. 

The layout of the paper is as follows. The simulations tech- 
niques and parameters are presented in Sect|2] distinguishing 
Tree-SPH and PM-SP models. Results from Tree-SPH simula- 
tions are presented in Sectl3] and, in Sect|4]the main findings are 



compared to that obtained employing a PM-SP approach. After 
comparing our results about the frequency of starburts episodes 
with observations (Sect|6ll, in Sect|7] the main conclusions of 
this research are presented. 

2. Description of the numerical simulations 

2.1. The data set 

In DM07, a statistical study of the relationship between star for- 
mation and galaxy interactions was performed using a set of 216 
simulations. The morphology of the interacting galaxies, as well 
as their encounter velocity and distance, were varied, and the role 
played by different parameters (gas fraction, galaxy minimum 
separation, galaxy relative velocity, strength of tidal effects, etc) 
was analyzed. That work represented also the first numerical 
work, to our knowledge, where the difficulty to drive substan- 
tial bursts of star formation during major galaxy encounters was 
outlined. The results presented in that paper showed, indeed, that 
major mergers were neither always responsible for triggering in- 
tense bursts of star formation, nor for converting large quantities 
of gas mass into new stars. However, the influence of the orbit 
inclination was not studied, all disk galaxies in this first sample 
being coplanar. Also, the influence of the numerical techniques 
and recipes to compute the star formation rate has not been stud- 
ied. The goal of the present work is to remove ffiese limitations. 

First, we extend the analysis of the relation between galaxy 
interactions and star formation to a larger parameter set, taking 
into account also encounters between galaxies having not zero 
disk inclinations. Indeed, coplanar interactions are quite pecu- 
liar cases that can favor the formation of regions of gas shocks, 
thus it is interesting to explore to what extent removing this as- 
sumption can modify the star formation evolution of the pair 
We will also study the star formation activity when the two in- 
teracting galaxies have an initial amount of gas higher than that 
of galaxies in the Local Universe. This initial work is here ex- 
tended with the realisation of a sample of 672 simulations (648 
for local galaxies and 24 for high-redshift systems) made with 
the same Tree- SPH code as in DM07, whose main characteris- 
tics are recalled below. 

To study whether the main conclusions depend on the numer- 
ical techniques and star formation schemes adopted, we have run 
a second set of 96 simulations with a different numerical code (a 
particle-mesh code wiffi a sticky-particle modeling of ffie ISM), 
where the initial conditions are also somewhat different. Some 
simulations in this sample also employ star formation models 
that differ from the Schmidt law. This is in order to understand 
how and in which way different star formation recipes can af- 
fect the star formation evolution during galaxy interactions, and 
if the conclusions of the main dataset are crucially dependent on 
the star formation model or not. 

In our opinion, the study of such a large and vast data set 
should help in finding strong and robust results about the star 
formation efficiency during galaxy major interactions, avoiding 
some of the limitations which have affected previous numerical 
works. 

2.2. Numerical methods 
2.2.1. Tree-SPH simulations 

The first set of 672 simulations has been realized employing ffie 
same Tree-SPH code as in DM07. Gravit ational forces are cal- 
culated using a hierarchical tree method dBarnes & HutI [19861) 
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and gas evolution is followed by means of smoothed particl e 
hydrodynamics (SPHj Lucvlll977l:lGingold & Monaghanlll982h . 
The code and the adopted numerical parameters are described 
in DM07 and references therein; we here recall only the main 
features. Gravitational forces are calculated using a tolerance 
parameter 6 - 0.7 and include terms up to the quadrupole or- 
der in the multiple expansion. A Plummer potential is used to 
soften gravity at small scales, with constant softening lengths 
£ - 280 pc for all species of particles. The equations of mo- 
tion are integrated using a leapfrog algorithm with a fixed time 
step of 0.5 Myr. A conventional form of the artificial viscosity 
is used in the SPH model, with parameters a - 0.5 and /? = 1 .0 
(iHern quist & Katz 1989). To describe different spatial dynami- 
cal ranges, SPH particles have individual smoothing lengths /i,, 
calculated in such a way that a constant number of neighbors lies 
within 2hi. The simulations in this paper have been performed 
using a number of neighbors A^, ~ 15. The interstellar gas is 
m odeled as isothermal , with a temperature Tga, = lO^/T. 
In lSemelin & CombesI (l2002h . a standard validation test for this 
code (the collapse of an initially static, isothermal sphere of self- 
gravitating gas) has been presented. No change is made in the 
algorithm with respect to that paper, and we here use an isother- 
mal equation of state for the gas instead of a multiphase model. 
Other tests on the dependency of the star formation rates on the 
numerical parameters adopted are presented in AppendixlAl 



2.2.2. PM-SP simulations 

The second set of simulations was performed with a PM-SP 
(particle-mesh - sticky- particles) code, which i s a gri d-based N- 
body code described in iBournaud & CombesI (l2003l) . The den- 
sity of particles is computed on Cartesian grids of maximal res- 
olution 330 pc through a Cloud-In C ell multi-linear interpola- 
tion. An FFT technique (|Jameslll977l) is then used to compute 
the associated gravitational potential, with a softening length 
e = 330pc. 

The collisional dynamics of the ISM is modeled with a 
sticky- particle scheme: gas particles, that model interstellar gas 
clouds, undergo inelastic collisions during which their relative 
velocity is multiplied by /3r along the direction of their positions, 
and yS, in the perpendicular direction. Here, we use Pr - 0.7 and 
jSf = 0.5, which creates a cool enough medium for thin spiral 
arms to form in isolated disks but without making the disk un- 
stable to axisymmetrical perturbations. In the case of particularly 
gas-rich galaxies (see Sect. 14. 3t , the dissipation is reduced, using 
jSr - 0.8 and f), - 0.7, in order to limit the formation of clump 
instabilities. 



2.3. Star formation recipes 

Different recipes and numerical methods have been adopted in 
existing models to include star formation and accoun t for the ef- 
fects that this star formation has on the surroundings (Katz' 1992'; 
Steinmetz & Muller 1994; Springel 2000; Springel & Hernquist 
20031; ICox et al.ll2006l) . These are usually based on the so-called 
Schmidt-Kennicutt law (Schmidt 1959, Kennicutt 1998a,b) and 
assume that the local star formation rate can be inferred from the 
local gas density, sometimes combined with a stability threshold. 
In the following two subsections, the diff'erent numerical recipes 
adopted in our codes are described. Their main features are sum- 
marized in Table [T| 



Table 1. Main features of the star formation recipes adopted in 
the Tree-SPH code and in the grid-sticky particles one. 







Tree-SPH 


PM-SP 




^SFR Sgas''^ 


yes 


yes 


SF laws 


^SFR '^gas^ 


no 


yes 






no 


yes 


SF feedback 


kinetic feedback 


yes 


no 




metal enrichment 


yes 


no 



2.3.1. Tree-SPH simulations 

In all the simulations performed adopting the Tree-SPH code, 
a density-dependent star formation law (Schmidt law) has been 
employed. The star formation rate is locally given by: 

PSFR^CXpgcJ-^ (1) 

with the constant C chosen such that the isolated disk galax- 
ies form stars at an average rate of between 1 and 2.5 M© yr"'. 
This parametrization is consistent with the observational evi- 
dence that on global scales the SFR in disk galaxies i s well rep- 
resented by a Schmidt law with an exponent 1 .4 (see Kennicud 
J998a b, but also Wong & Blitz 2002; Boissier et al. 20031 
iGao & Solomonll2004 . 

Once the SFR recipe is defined, we ap ply it to SPH parti- 
cles, us ing the "hybrid" method described in lMihos & Hernquist! 
(Il994bh . This consists in representing each gas particle with two 
mass values, one corresponding to its gravitational mass M,, 
whose value stays unchanged during the whole simulation, and 
the other describing the gas content of the particle M, ^aj, whose 
value changes over time, according to EqUI Gravitational forces 
are always evaluated on the gravitational mass M,, while hydro- 
dynamical quantities, in turn, uses the time-varying mass of gas 
Mi^gas- If the gas fraction present in the hybrid particles drops 
below 5% of the initial gas content, the hybrid particle is totally 
converted into a star-like particle and the small amount of gas 
material still present is spread over neighbors. 

We also fo ll owed the method described in 
iMihos & Hernquisj (Il994bl) for including the effects of 
star formation into the interstellar medium (metal enrichment 
and energy injection in the ISM by supernovae explosions). 
The method is fully described by these authors, so we refer the 
reader to it for more details. 

2.3.2. PM-SP simulations 

The Tree-SPH simulations have been carried out for a large set 
of morphological and orbital parameters (see Sect. 12.4. It , but 
only with a Schmidt law to compute the star formation rate. The 
set of orbital parameter in the PM-SP sample is less extensive 
(Sect. |2.4.2t , which enabled us to perform simulations with dif- 
ferent models to compute the star formation rate: 

- First, we used the same Schmidt law as in Tree-SPH simu- 
lations with an exponent 1.5, as described in Eq.[r| There is, 
however, a diff'erence in the way this prescription is applied. 
In Tree-SPH simulations, the gas density pgas is the SPH den- 
sity, computed with an adaptative resolution. Some theory 
predicts that the Schmidt law is scale-free (e.g., Elmegreegl 
2002) but this is not necessarily the case. In the PM-SP sim- 
ulations, we then compute the density pgas on the Cartesian 
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Fig. 1. Comparison of the three star formation models used in the PM-SP simulations to observational data from lKennicuttI ( Il998bh . 
Left: Schmidt law with a 1 .5 exponent - Middle: 'Ls fr '^gas x ^ model, also compatible with observations - Right: l,s fr '^I'as ^ ^ 
model. This last model is barely compatible with observations and tends to overestimate the star formation rate at high gas density, 
we thus use this model to put an upper limit on the actual starburst efficiencies. 



Table 2. Tree-SPH and PM-SP simulations: Galaxy parameters. The bulge and the halo are modeled as Plummer spheres, with 
characteristic masses Mb and Mh and characteristic radii and rn- M, and Mg represent the masses of the stellar and gaseous 
disks, whose vertical and radial scale lengths are given, respectively, by and and lig and Og. 





gEO 


Tree-SPH 
gSa gSb gSd 


gSb+ 


gSbO 


PM- SP 

gSb 


gSh+ 


Mb [2.3 X 10" Mo] 


70 


10 


5 





5 


5 


5 


5 


Mh [2.3 X 10" Mq] 


30 


50 


75 


75 


75 


65 


65 


65 


M. [2.3 X 10" Mo] 





40 


20 


25 


20 


22 


20 


12 


Mg/M, 





0.1 


0.2 


0.3 


0.5 





0.15 


0.50 


rg [kpc] 


4 


2 


1 




1 


1.8 


1.8 


1.8 


rH [kpc] 


7 


10 


12 


15 


12 


10. 


10. 


10. 


a, [kpc] 




4 


5 


6 


5 


5. 


5. 


5. 


h, [kpc] 




0.5 


0.5 


0.5 


0.5 


0.7 


0.7 


0.7 


Og [kpc] 




5 


6 


7 


6 


15. 


15. 


15. 


hg [kpc] 




0.2 


0.2 


0.2 


0.2 


0.25 


0.25 


0.25 


Qgas 




0.8 


0.8 


0.8 


1. 


1.2 


1.2 


1.2 



Table 3. Tree-SPH and PM-SP simulations: Particle numbers for each galactic component 







Tree-SPH 




PM- 


SP 




gEO 


gSa gSb, gSb+ 


gSd 


gSbO i 


;Sb,gSb+ 






20000 40000 


60000 




100000 




80000 


60000 40000 


20000 


70000 


70000 


Ndm 


40000 


40000 40000 


40000 


50000 


50000 



grid with a fixed resolution. The corresponding ps fr indi- 
cates a number of gas particles to be converted into stellar 
particles in each cell. 
- Second, we model a star formation rate psfr °^ pgas x 
where Q. is the local rotational angular velocity of the gas 
disk. This can model the smaller size and short er collapse 
timescale of molecular clouds at small radii (Elmegreen 
ri997) and/or their higher collision frequency (.SiUc 1997). 
Just like the Sch midt law, th i s is co mpatible with observa- 
tions according to Kennicutl ( 1199 8b') - see also Fig. [T] It is 
here practical that only one galaxy contains gas in the PM-SP 



simulations (see Sect l2.4.2b . so that the time-dependent rota- 
tion curve Q(r) can be computed for this galaxy. To this aim, 
we compute the velocity of gas particles in 50 pc wide ra- 
dial bins every 75 Myr (assuming Q. is constant in the central 
100 pc to avoid singularities). We use the result as a proxy for 
the actual rotation curve, which could anyway not be clearly 
defined during the merger This is a simple way to estimate 
the influence of the dynamical timescale on the star forma- 
tion activity during the merger. We checked that using the 
initial rotation curve of the galaxy instead does not lead to 
large changes: the rotation curve evolves mainly during the 
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Fig. 2. Adopted orbital geometry for our simulations. We set up 
the collision in such a way that the orbital angular momentum 
is parallel to the z-axis and that the centers of the two galaxies 
are on the x-axis initially. The galaxy spins are represented by 
the blue and red arrows, respectively. They are specified in terms 
of the spherical coordinates (i'l, <l>i) and (12, 02). See Table|4]for 
their initial values. 

very final stages of the merger relaxation, while most star 
formation occurs earlier, so that the way Q(r) is estimated is 
not crucial. 

- Third, we used a model with ps fr Pgas ^ which is the 
combination of the two previous models. This scheme has 
no theoretical motivation but corresponds to an upper limit of 
the observed non-linearity of the gas density - star formation 
rate density relation (see Figure [T]!. This model should thus 
correspond to an upper limit to the star formation efficiency 
in galaxy mergers. 

These three star formation recipes are applied in 3-D mod- 
els, but are equivalent to 2-D laws based on the gas surface den- 
sity E and SFR surface density, under the reasonable assump- 
tion of a uniform thickness of gas disks. Ea ch star formation 
model is compared to observational data from iKennicuttI (Il998bh 
in Fig.IU 

In the PM-SP simulations, we do not use hybrid particles like 
in the Tree-SPH model, but instead convert gas particles into 
stellar particles at a rate that, in each cell of the grid, is given 
by the chosen star formation model. Not using hybrid particles 
reduces the mass resolution in the treatment of star formation, 
but prevents a newly formed stellar mass to follow an SPH dy- 
namics. Each method thus has drawbacks, and the two different 
assumptions are tested in our two datasets. 

Energy feedback from supernovae is not included in the PM- 
SP simulations. The reasons for this choice are: (i) this provides 
different assumptions compared to the Tree-SPH sample that in- 
cludes feedback; ( ii) within the most standard assumptions, feed- 
back does not ha ve a major influe nce on the relative SFR evolu- 
tion in mergers jCox et alj|2006b and (Hi) the general effect is 
to regulate star formation and reduce the starburst efficiency in 
the merging phase- because our conclusion is that most merger- 
driven starbursts have a low efficiency, not including feedback is 
a rather conservative choice. 

2.4. Initial conditions 
2.4.1. Tree-SPH simulations 

Galaxy models: moving along the Hubble sequence 

As in DM07, our aim is to exploit a large set of interactions, 
involving galaxies of all morphologies from ellipticals to late- 
type spirals. For each galaxy type, the halo and the bulge (if 



Table 4. Orientation of the galaxy spins, for Tree-SPH and PM- 
SP simulations. 





Tree-SPH 


PM-SP 




0° 


33° 


0, 


0° 


30° 




0°, 45°, 75°, 90° 


0° 


O2 


0° 


0° 



any) are modeled as a Plummer sphere jBinnev & Tremain3 
11987', pag.42), with characteristic masses Mb and Mh and 
characteristic radii rg and rn. The stella r and gaseous disks 
follow a Miyamoto-Nagai density profile jBinnev & Tremaing 
1987, pag.44) with masses M, and Mg and vertical and radial 
scale lengths given, respectively, by /i* and a*, and hg and ag. 
For the different morphologies adopted, stellar disks have a 
peak rotation speed within and 3fl, which is between 63% 
(for Sbc spirals) and 75% (for Sa and Sd spirals) of the total 
circular speed at that radii. The initial rotation curves of these 
models are given in DM07, Fig. 2. 

Hereafter, we will adopt the following nomenclature for the 
different morphological types: gEO for giant-like ellipticals, gSa 
for giant-like Sa spirals, gSb for giant-like Sbc spirals and gSd 
for giant-like Sd spiral£l For giant-like Sbc spirals, we will also 
perform some simulations adopting a higher gas mass fraction 
than that typical of galaxies in the Local Universe. We will use 
the nomenclature gSb-i- to distinguish these galaxies with a high 
gas mass fraction from the "local" Sbc ones. The complete list 
of all the parameters is given in Table |2] We refer the reader to 
DM07, for a representation of our galaxy sequence. 

Since we investigate interactions between giant-like galax- 
ies, the mass ratios of the interacting systems is always of the 
order of unity. Intending to run hundreds of simulations, each 
galaxy is made up of 120000 particles, distributed among gas, 
stars and dark matter, depending on the morphological type (see 
Table|3]l. To i nitialize particle velocities, we adopted the method 
described in iHernquislI (119931) . 

Orbital parameters 

In DM07, seeking to exploit a vast range of orbital parameters, 
we performed 24 different simulations for each couple of inter- 
acting galaxies, varying the orbital initial conditions, in order to 
have (for the ideal Keplerian orbit of two equal point masses of 
mass m - 2.3 x 10"Mo) the first pericenter separation r^g,. = 
8, 16, and 24 kpc. For each of these separations, we varied the 
relative velocities at pericenter, in order to have one parabolic 
and three hyperbolic orbits of different energy. Finally, for each 
of the selected orbits, we changed the sign of the orbital angu- 
lar momentum in order to study both direct and retrograde en- 
counters. Combining each orbital configuration with all possible 
morphologies for the interacting pair of galaxies, we obtained a 
total sample of 216 interactions, including only coplanar pairs. 

In this work, we have extended the previous sample, taking 
into account also different disk inclinations for the interacting 
galaxies. In particular, for each interacting pair in the Tree-SPH 
sample, we have kept the disk (when present) of one of the galax- 
ies in the orbital plane (/] - 0°), and varied the inclination (2 



' We use the term "giant" for these systems, to indicate that they are 
not dwarf galaxies. Their masses are indeed comparable to that of the 
Milky Way. 
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Table 5. Tree-SPH simulations: Galaxies orbital parameters 



id 


^ ini 


^/ 






spin"^ 




[kpc] 


[kpc] 


[lO^kms-'] 






Oldir 


100. 


8.0 


7.07 


0.0 


up 


Olret 


100. 


8.0 


7.07 


0.0 


down 


02dir 


100. 


8.0 


7.42 


2.5 


up 


02ret 


100. 


8.0 


7.42 


2.5 


down 


03dir 


100. 


8.0 


7.74 


5.0 


up 


03ret 


100. 


8.0 


7.74 


5.0 


down 


04dir 


100. 


8.0 


8.94 


15.0 


up 


04ret 


100. 


8.0 


8.94 


15.0 


down 


05dir 


100. 


16.0 


5.00 


0.0 


up 


05ret 


100. 


16.0 


5.00 


0.0 


down 


06dir 


100. 


16.0 


5.48 


2.5 


up 


06ret 


100. 


16.0 


5.48 


2.5 


down 


07dir 


100. 


16.0 


5.92 


5.0 


up 


07ret 


100. 


16.0 


5.92 


5.0 


down 


08dir 


100. 


16.0 


7.42 


15.0 


up 


08ret 


100. 


16.0 


7.42 


15.0 


down 


09dir 


100. 


24.0 


4.08 


0.0 


up 


09ret 


100. 


24.0 


4.08 


0.0 


down 


lOdir 


100. 


24.0 


4.65 


2.5 


up 


lOret 


100. 


24.0 


4.65 


2.5 


down 


lldir 


100. 


24.0 


5.16 


5.0 


up 


llret 


100. 


24.0 


5.16 


5.0 


down 


12dir 


100. 


24.0 


6.83 


15.0 


up 


12ret 


100. 


24.0 


6.83 


15.0 


down 



" For two equal point masses with in = 2.3 x 10" Mq. 

It is the total energy of the relative motion, i.e. 

E = - G(mi + m2)/r. 
° Orbital spin, if parallel (up) or antiparallel (down) to the galaxies 

spin. 



of the companion disk, considering : 12 - 45°, 12 - 75° and 
12 - 90°, for a total of 648 new simulations (see Fig|2] for a 
sketch of the initial orbital geometry and Table |4] for the coordi- 
nates of galaxy spins). 

In Table |5] the initial distance r„„ and the pericenter distance 
rp between the galaxies center-of-mass are listed, together with 
their relative velocity Vp at pericenter and the orbital energy E, 
for all the simulated encounterfl 

We have also run 24 additional simulations of two interacting 
coplanar (/i = 12 - 0°) gSb+ galaxies, for all the possible orbital 
parameters given in Table |5] 

In the following, we will refer sometimes to specific encoun- 
ters, by adopting the morphological type of the two galaxies in 
the interaction (gEO, gSa, gSb or gSd), + the encounter identifi- 
cation string (see first column in Table [S]), + disk inclination 12 
of the second galaxy (ii is always equal to zero). For example, 
the nomenclature gSagSb04ret45 coiTesponds to an interaction 
between a giant Sa and a giant Sb spiral; the Sb disk is inclined 
of i2 - 45 degrees with respect to the orbital plan^ and the ini- 
tial orbital parameters of the encounter are those corresponding 
to id=04ret in Table |5l 

2.4.2. PM-SP simulations 



^ The values refer to the ideal Keplerian orbit of two equal point 
masses of mass in = 2.3 x 10"Mg. 

For Tree-SPH simulations, ii being always null, the relative incli- 
nation between the two galaxy disks is equal to !2. 
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Fig. 3. Face-on and edge-on views (left : stars, right : gas) of the 
"evolved" initial conditions in the PM-SP simulations, here for 
the gSb model galaxy. Each box is 30 kpc x 30 kpc in size. 

Galaxy models 

The galaxy mergers in DM07 and in the larger Tree-SPH sam- 
ple described above do not show a drastic evolution of the star 
formation efficiency with the internal properties of the merging 
galaxies. In the sample of PM-SP simulations, we model only 
gSbO-gSb encounters, which are basically representative on av- 
erage of any galactic encounter 

The gSb-PM model galaxy is roughly similar to the gSb model 
from the Tree-SPH sample. The stellar disk has a Toomre ra- 
dial profile, the bulge and dark matter halo a Plummer profile, 
with mass and scale-lengths indicated in Table|2] A radial scale- 
length for gas larger than that of Tree-SPH simulations is used. 

The gSbO-PM model has the same parameters than the gSb- 
PM galaxy, except for the fact that the disk does not contain gas. 
In other words, all the PM-SP runs concern the interaction of 
two spiral galaxies, with a bulge-to-disk ratio equal to 0.2, only 
one of the two disks containing gas. 

The galaxy interactions in the Tree-SPH simulations are 
started with axisymmetric disks (see, for example, Fig.7 in 
DM07). This is common practice in such numerical studies, 
however, the rapid transition from these initial conditions to a 
more realistic (frequently baiTed) spiral distribution occurs dur- 
ing the early phases of the interaction/merger. To make sure that 
the future conclusions are really related to the merger-driven 
evolution and not to this artificial evolution of the initial con- 
ditions, a different method was used in the PM-SP simulations. 
The simulated galaxies have been evolved during 1 Gyr in iso- 
lated conditions, before the simulations of interaction are started. 
In this way, the interacting galaxies already have a realistic and 
slowly evolving barred spiral structure in the distribution of their 
stellar and gaseous components: this helps to make sure that the 
SFR evolution relates to merger-driven processes and not to the 
rapid formation of the initial spiral structure. This initial 1 Gyr 
evolution is made over a short period compared to secular evolu- 
tion timescales, so that the bulge mass or disk scale-length do not 
change dramatically compared to their initial values. However, 
star formation is turned off during this initial phase, so that the 
gas fraction we have indicated (15 %) really corresponds to the 
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gas fraction in the galaxies when the interaction starts, and any 
spurious consumption of gas during the initial formation of the 
bar and spiral arms is avoidecfl We show these "evolved" initial 
conditions of the PM-SP models on Fig. [3] 



Orbital parameters 

The analysis of the Tree-SPH sample (below) shows that the 
inclination does not have a major impact on the statistics of star- 
burst duration and intensity. We thus choose not to vary this pa- 
rameter in the PM-SP sample, in order to reduce the number of 
simulations and leave the possibility to vary the star formation 
scheme instead. The inclinations are then fixed to ii = 33° and 
12 = 0° , because this is about the average value for an isotropic 
distribution of mergers (Bournaud et al. 2007h) . The other or- 
bital parameters are varied as follows: 

- the encounter velocity V to 50, 100, 150 and 200 km s"' 

- the impact parameter b to 37.5, 50, 62.5 and 75 kpc 

- the orientation to prograde and retrograde. 

Each combination of these parameters that eventually leads to a 
merger of the galaxy pair has been simulated and is included in 
the following statistical analysis. 

2.4.3. Core versus cuspy dark matter profiles 

In both Tree-SPH and PM-SP sets of simulations, the dark mat- 
ter profiles have been modeled with a core density distribution, 
which seems to be more in accordance with the rotation curves 
of local spirals and dwarf galaxies rather than the cuspy halos 
predicted by Cold Dark Matter (CDM) numerical simulations. 
Indeed, while cosmological numerical simulations of CDM pre- 
dict cuspy profiles with a density distribution showing a p ~ r'"' 
behav i or, wi t h a > 1, (jC ole & Lacey 1996; Navar ro et alj 



2005D, 



119961 119971: lAvila-Reese et all 119981: iDiemand et all 
numerous studies of the rotation curves of local galaxies 
have shown the observational data are more consistent with 
the presence of a dark halo having a nearly flat density 
core (a ~ 0). T his has been shown to b e valid either for 
dwarf galaxies (Flores & PrimacS 119941 : iMarchesini et alj 
I2OO2I: iGentile et al. .2005). or for L ow Surface Brightness 



(LSB) galax ies dde Blok et al J 1200 U iM cGaugh et al. 2001; 
IMarchesini et alj 



200a iKuzio de Narav et al., i2006t i2008l) . 



Some authors have driven the attention to the possibility that 
systematic effects in the data (as beam smearing, non-circular 
motions, slit inclina tions) c ould mask the presence of a cusp 
dvan den B osch all I2OO0I: Ivan den Bosch & S waters! 1200 U 
ISwaters et alj i2003l) . suggesting a possible bimodality in the 
dark halo profiles, with constant density cores restricted to 
low mass system s (Vmax < lOkmJ s) a nd cuspy profiles for 
hi gh mass ones (ivan den Bosch aljr2000l) . However, according 
to ldeBloket~an (l2003l) none of these systematic effects can 
reconcile the observational data with cuspy CDM halos. In any 
case, a s shown for the LSB galaxy UGC4325 by Che min et alj 
(I2007h . an accurate modelisation of gas motions seems to be 
necessary in order to disentangle the dark halo shape of these 
systems. 



In the case of Tree-SPH runs, in turn, simulations are started with 
axisymmetric disks and star formation is included since the beginning. 
Thus, in the merging phase, galaxies have yet consumed a fraction of 
their gas content (typically 50 %, even if with large variations, as shown 
in Fig. 20 of DM07). Nevertheless, as discussed in that paper, this does 
not lead to systematic effects on the relative star formation rate of the 
pairs. 



Moreover, recent re sults from a new HI mapping of M31 
dChemin et alj 12008 ^ also seem to rule out the presence of a 
Navarro-Frenk- White cusp in the dark halo profile of this galaxy. 

The adoption of a cuspy dark halo profile and the induced 
star formation in merging systems will be the subject of a fu- 
ture paper (Combes et al, in preparation); in any case, some 
comparisons of bars and gas evolution in isolated galaxies with 
cuspy Navarro-Frenk- White dark halo profiles or core Plummer 
ones are discussed in ICombesI (l2007h . These simulations show 
that the dark matter profile plays an important role in the de- 
velopment of bar instabilities, and that the presence of a cuspy 
profile tends to create axisymmetric mass concentrations, which 
dilutes the gravity torques of the bar and the subsequent gas 
flows. Similarly, DM cusps may reduce the merger-driven in- 
flows (the driving mechanism, gravity torques, being the same 
as for bar-driven flows), which in turn may somewhat decrease 
the resulting star- formation rate. Our choice of core haloes is 
then in better agreement with observations of local spirals, and is 
also a conservative hypothesis regarding the intensity of merger- 
induced starbursts. 



3. Results from the Tree-SPH simulations 

3.1. Some SFR evolutions 

In this section, we present some SFR evolution during galaxy 
interactions. In order to distinguish the contribution of tidal ef- 
fects from secular evolution in determining the star formation 
history of the two galaxies, all the star formation rates presented 
in this and in the next sections are always computed relative to 
that of the two corresponding isolated galaxies. For the SFRs 
of the isolated galaxies, they have been discussed and shown in 
Sect.4.2 and Fig. 6 of DM07. In the case of Tree-SPH models, 
galaxy simulations are started with axisymmetric disks, so that 
the evolution of the isolated systems shows a transient initial 
burst which is due to compression of the gas into density waves. 
This transient burst, in turn, is not present in PM-SP simulations: 
in this case simulated galaxies have been run in isolated condi- 
tions for about 1 Gyr, before the interactions are started. 
FigHjshows some SFR evolutions with time during major merg- 
ers involving galaxies of different morphologies and different or- 
bital parameters. In general, there is a tendency to an enhance- 
ment in the star formation rate during an interaction, but, as it 
can be seen, a variety of star formation histories are found. In 
some cases, as in the encounter gSagSa05ret45 (first panel of 
Fig.©, the enhancement of the star formation rate is practically 
null at the first pericenter passage, while the SFR in the merg- 
ing phase peaks at more than 10 times that of the corresponding 
isolated case. The encounter of two gSa galaxies, one having 
a disk inclined of 90° with respect to the orbital plane (second 
panel of Fig. |4]l, leads to two peaks of star formation: the first, 
at the pericenter passage, being more intense than the one which 
takes place in the coalescence phase of the two galaxies. The 
interaction of a gSa and a gSb galaxies, on retrograde orbits, 
with a relative disk inclination of 45° (displayed in Fig.|4] third 
panel) causes, in turn, a SFR enhancement which starts at the 
first pericenter passage and lasts until the final merging phase. 
Even more striking is the case shown in the last panel of Fig.|4] 
concerning a giant Sa and a giant Sd moving on a direct orbit 
(gSagSd05dir45): in this case, the maximum in the star forma- 
tion evolution takes place neither at the pericenter passage, nor in 
the merging phase, but about 1 .6 Gyr after the coalesc ence of the 
two systems. At that epoch, a tidal dwarf galaxy (see lDuc et alj 
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Fig. 4. Star formation rate, versus time, for some galaxy mergers. The SFR is normaHzed to that of the corresponding isolated 
galaxies. In each panel, the black arrows indicate, respectively, the first pericenter passage between the two galaxies and the merger 
epoch. 
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Table 6. Number of mergers and flybys in 3 Gyr, the simulated 
time interval, for the different inclinations /. 



Fig. 5. Upper panels: gas+new stars maps for the remnant of the 
merger gSagSd05dir45. Lower panels: old stars maps. From left 
to right, maps are shown at t=1500 Myr, t=1550 Myr, t=1600 
Myr, t=1650 Myr and t=1700 Myr after the coalescence of the 
gSa and gSb galaxies. Note the presence of a dwarf galaxy, 
formed during the galaxy encounter, which is falling into the 
central region of the remnant. Each box size is 12 kpc in length. 



|2004|) . formed during the first phases of the interaction, falls into 
the central regions of the remnant, stimulating another episode 
of star formation (see Fig|5]for some maps of this remnant, dur- 
ing the accretion of the dwarf satellite). 

3.2. Starburst frequency 

In the previous section, we have shown that galaxy encounters 
can lead to a variety of star formation evolutions: some systems 
show SFR typical of starbust galaxies, while others show only 
a weak enhancement in the star formation rate during the inter- 
action. In this and in the following sections, we want to deepen 
the discussion on this point, because of its potential impact on 
observational and cosmological studies. 

In particular, in this section, we present some histograms of the 
maximum star formation rate for mergers and flyby^ As done 
previously, the SFR is relative to that of the corresponding galax- 
ies evolving isolated. Fig|6]shows these histograms for three dif- 



' By mergers, we mean an encounter which leads to the coalescence 
of the two galaxies in the simulated time interval (i.e. 3 Gyr), while a 
flyby is an encounter which does not lead to the coalescence of the two 
systems during the same time interval. 
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ferent groups of encounters: those having inclination of 12 = 45°, 
12 = 75° and 12 = 90° between the disk of one of the two galaxies 
and the orbital plane. For comparison, coplanar encounters are 
also shown (the results for coplanar encounters are extensively 
discussed in DM07). 

It results that mergers do not always trigger starbursts. 
Indeed the fraction of merging galaxies that produce star forma- 
tion rates at least ten times higher than those of isolated galaxies 
is about 17% of the total merger sample for encounters with in- 
clination 12 - 0°, 22% for encounters with inclination 12 =45°, 
while this fraction decreases to 15% and 13% for encounters 
with inclination of 12 - 75° and 12 - 90°, respectively. Moreover, 
as the inclination increases, the maximum amplitude for the star 
formation in the merging phase decreases: for coplanar merg- 
ers, the maximum SFR was about 60 times those of the isolated 
galaxies, while, for polar encounters, it is only 30 times greater 
than the star formation rate of the isolated systems. Note also 
that in this figure, as in the case of Figs|2l[8]and|9] the number of 
mergers Nmer and flybys Njiy changes with the inclination, while 
the total number of encounters Ntot - Nmer + Nfiy is equal to 
216, for all the inclinations chosen (see Table|6]l. 

3.3. Duration of the star formation enhancement 

Intense starbursts during galaxy mergers are not only less fre- 
quent, but also characterized by shorter duration times. In Fig|7] 
we present the duration of the star formation enhancement dur- 
ing galaxy mergers and flybys. More precisely, we investigate 
the duration of the star formation enhancement above a certain 
threshold, which has been taken equal to two and five times the 
star formation rate of the isolated galaxies. The histograms in 
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Fig. 6. Histograms of the maximum SFR (relative to the isolated case) for mergers. Flybys are shown for comparison in the small 
window inserted in the figure. From left to right: histograms relative to encounters with 12 - 0°, 12 = 45°, 12 = 75° and 12 - 90°, 
respectively. For each inclination, the total number of mergers and flybys is shown in Table |6] More details on the statistical 
distribution of the maximum SFR are given in Table IdTT] AppendixiDl 
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Fig. 7. Histograms of duration of enhanced SFR for the whole sample of interacting galaxies (mergers and flybys). Two thresholds 
are shown: relative SFR > 2 and relative SFR > 5 (the latter is shown in the small window inserted in the figure). From left to 
right: histograms relative to encounters with 12 - 0°, 12 - 45°, 12 - 75° and 12 - 90°, respectively. More details on the statistical 
distribution of the duration of enhanced SFR are given in Table ID. 2l AppendixiDl 
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Fig. 8. Histogram of the relative ISFR for mergers, Flybys are shown for comparison in the small window inserted in the figure. 
From left to right: histograms relative to encounters with 12 - 0°, 12 = 45°, 12 = 75° and 12 = 90°, respectively. For each inclination, 
the total number of mergers and flybys is shown in Table |6] More details on the statistical distribution of the ISFR are given in 
Table |D31 Appendix |D] 



Fig|2]refer only to encounters which sustain a star formation en- 
hancement greater than the chosen thresholds, 

- Among the 216 encounter^ with inclination 12 - 45°, 185 
sustain a star formation rate which is at least 2 times higher 
than that of their isolated counterparts (the remaining 3 1 en- 
counters show no significant enhancement). Among these 
185 encounters, only 26 (i,e, only 14%) are able to sustain 
this star formation rate for more than 500 Myr, 

- Among the 216 interactions with 12 =45°, only 46 show a 
relative star formation enhancement greater than 5 (i,e, about 
21% of the total sample). Among these 46 encounters, 20 of 
them sustain this enhancement for a time greater than 100 
Myr 

* i,e, including mergers and flybys. 



In other words, according to these results, in an interacting 
galaxy sample, the probability of finding starburst galaxies rather 
than "normal" ones should be small, not only because starbursts 
are less frequent, as we saw in Sect l3,2l but also because the 
duration of the star formation enhancement shortens as the rate 
of star formation increases. Similar trends are found also in the 
two panels in Fig|7] which refer to encounters with inclination 
i2 = 75° andi2 =90°, 

3.4. Integrated star formation rate 

The results presented in the previous sections suggest that the 
relation between galaxy interactions and star formation is quite 
complex. Some systems show in the merging phase a star for- 
mation rate which can be ten times higher than that of isolated 
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Fig. 9. Histogram of the relative ISFR for mergers. Flybys are shown for comparison in the small window inserted in the figure. 
From left to right: histograms relative to encounters with 12 = 0°, 12 = 45°, 12 = 75° and 12 = 90°, respectively. The top row refers 
to the ISFR between t\ = tper - 50 Myr and t2 = tp^r + 300 Myr, being tp^r the time of the first pericenter passage. The bottom row 
shows the ISFR in a time f between f2 - tper + 300 Myr and 3 Gyr. For each inclination, the total number of mergers and flybys is 
shown in Table |6|. 



galaxies, while many systems show only a weak enhancement. 
We have also seen that the duration of the enhanced phase of star 
formation crucially depends on the level of SFR sustained: the 
higher the star formation rate, the lower its duration. The next 
step is to quantify the total gas mass converted into stars dur- 
ing an interaction. As in DM07, we quantify the integrated star 
formation rate during a time interval T - (t\, t2) as 



ISFR 



nt=t2 nt=t2 

I SFR{t)dtl I SFRiso(t)dt 

Ji=i, Jt=t, 



(2) 



SFR(t) being the star formation rate of the interacting pair 
at time t and SFRi^o that of the two corresponding galaxies 
evolving isolated. In this way, it is possible to distinguish secular 
evolution from the effects due to the tidal encounter 
When integrating over the whole duration of the simulations 
(i.e. over 3 Gyr), independently of the disk inclination (2, we 
find that mergers do not always convert high gas mass quantities 
into new stars (see Figl8]l. While coplanar encounters can be 
efficient enough to produce twice as many stars as isolated 
galaxies, varying the disk inclination decreases the maximum 
amount of gas mass transformed into stars: in all cases, merging 
galaxies can produce, in the most favorable cases, 1.7 times 
more stars than their isolated counterparts. But the bulk of 
merging galaxies shows only a modest enhancement in the 
ISFR: for i2 = 45°, for example, 74% of the total merger sample 
show an ISFR < 1.3 times the integrated star formation of 
isolated galaxies. It is also interesting to note that, as the disk 
inclination increases, moving from ;2 - 45° to i2 = 90°, the 
number of mergers which shows a high value of the ISFR 
decreases. For example, mergers with an IS FR > 1.5 times the 
IS FR of the isolated galaxies constitute 14% of the total merger 
sample with i2 = 45°, while they are about 7% for mergers with 
!2 = 75° and / = 90°. 



In order to understand when most of the gas is transformed 
into stars during an interaction, we also analyzed the ISFR for 
two different time intervals: 

1 . the first includes the first phase of the interaction, going from 



ti 



^per 



-50Myr to f2 



'■per 



- 300Myr, where tp^r is the time 



of the first pericenter passage between the two galaxies; 
2. the second time interval includes the subsequent phase of the 
interaction, going from t\ - tp^r + 300Myr to t2 - 3Gyr, the 
end of the simulation. 

We want to emphasize that, while for mergers, the second 
time interval includes the coalescence phase of the two galaxies, 
for flybys most of the tidal effects act when the galaxies 
approach one another (so in the first time interval). 
The results of this study are shown in Fig|9] When looking at 
mergers, not surprisingly, we find that most of the conversion 
from gas to new stars takes place in the second phase of the 
interaction, for tper + 300Myr < f < 3 Gyr. Anyway, it should be 
noted that most of the galaxies in the late phases of the merging 
process (i.e. in the second time interval chosen) shows only a 
moderate enhancement in the ISFR: for inclination i2 - 45°, for 
example, 69% of the merger sample show an ISFR which is less 
than 1.3 times that of the corresponding isolated galaxies, and 
only 19% increase their ISFR of a factor greater than 1.5 with 
respect to isolated systems. 



Finally, it is interesting to note the response of galaxies to 
high velocity encounters, which do not lead to the coalescence 
of the two systems (flybys). Quite surprisingly, these encounters 
show an enhanced ISFR also in the second time interval, i.e. well 
after the pericenter passage, when the systems are at this point 
well separated. This means that galaxy encounters can stimulate 
an increase in the star formation rate of a galaxy not only during 
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Fig. 10. Gas maps for the isolated gSb+ galaxy. From left to 
right, maps are shown from t=500 Myr to t=3 Gyr, every 500 
Myr Both xy projection (top panels) and xz projection (bottom 
panels) are shown. Each box is 40 kpc x 40 kpc in size. 
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Fig. 11. Star formation rate, versus time, for some coplanar 
mergers involving two gSb+ galaxies. The SFR is normalized 
to that of the corresponding isolated galaxies. In each panel, the 
black arrows indicate, respectively, the first pericenter passage 
between the two galaxies and the merger epoch. 
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Fig. 12. Star formation rate, versus time, for some coplanar 
mergers involving two gSb galaxies. The SFR is normalized to 
that of the corresponding isolated galaxies. In each panel, the 
black arrows indicate, respectively, the first pericenter passage 
between the two galaxies and the merger epoch. See FigfTTjfor 
a comparison of these SFRs with those of gSb+ galaxies, having 
the same orbital parameters. 

more unstable, and the instability of their di sks can give rise 
to int ense star formation (see for example lElmegreen et al.l 
( |2QQ7|)), via fragmentat ion of the gas disk and clumps formation 
(iBournaud et al.ll200 7a). We also kept unchanged the extension 
of the stellar and gaseous disks, even if observations suggest that 
the si z es of disks decrease sign i ficantly with redshift ( Mao et aD 
I1998L iGiallongo et aP Il999t lAvila-Reese & Firmanii i200li 
These choices are all consistent with the fact that we want to 
investigate the effect of increasing the gas fraction in the disk, 
rather than accurately modeling galaxies at higher redshifts. 



the phases of close passage, but also well after this phase, when 
the two systems are far away from each other This delayed star 
formation can be caused by the infall of satellite dwarf galaxies 
formed during the interaction, as shown in FigH) or from insta- 
bilities in the galaxy disks, caused by the encounter. 

3.5. Increasing the initial gas fraction 

All the results presented in the previous sections concern the 
interaction of galaxies having a gas mass content typical of 
galaxies in the Local Universe. In this section, we want to 
discuss how and to what extent our results depend on the total 
amount of gas available in the galaxy disk. In other words, how 
does the star formation rate change, when the gas mass fraction 
in the galaxy increases? To answer this question, we have run 
a small number of simulations (24) of coplanar interactions 
between two Sb galaxies, having an initial gas mass equal to 
one half that of the stellar dislQ- 

The initial conditions chosen for the gSb+ galaxies are such 
that the initial disks are stable (see Fig[TO] for some gas maps 
of the gSb-H galaxy, evolving isolated for 3 Gyr). For that, 
we have increased the initial value of the Toomre parameter 
Qgas (see Table Probably, galaxies at higher redshifts are 



' Hereafter called gSb+ galaxies, to distinguish them from the gSb 
galaxies, presented so far. 



Some evolutions of the star formation rate, as a function 
of time, for gSb+-gSb+ interactions, are given in Fig[TT] As 
previously done, the SFRs are normalized to those of the 
corresponding isolated galaxies (but see Appendix |B] for the 
evolution of some absolute SFRs). 

The left panel in this figure presents the star formation evolution 
during a retrograde merger (id=gSb+gSb+01retOO): in this case, 
the coalescence of the two galaxies takes place only 200 Myr 
after the first pericenter passage. The stimulated star formation 
rate, after the first encounter, is quite high (about 15 times that 
of the isolated galaxy), while its amplitude is reduced in the 
merging phase, when the SFR is only 4 times that of the isolated 
counterparts. 

The response of the star formation evolution to tidal effects 
is also shown in the case of the retrograde merger with 
id=gSb+gSb+09ret00 (right panel in FigfTTTl. In this case, after 
the first passage (t=400 Myr), the star formation rate increases 
up to about 3.5 times that of the corresponding isolated galaxies 
(t=600 Myr), then it declines to preinteraction levels, to rise 
finally in the merging phase up to about 3 times the correspond- 
ing value of the isolated systems. 

Comparing these SFR evolutions with those of gSb galaxies, 
having the same initial orbital conditions, quite surprisingly we 
find that in the merging phase the SFR (relative to the isolated 
case) is higher for local gSb galaxies, than for gSb-H systems, 
which have initially a higher gas mass fraction (cf, for example. 
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Fig. 13. Gas maps (top) and star forming regions (bottom) for a retrograde merger involving two gSb+ galaxies 
(id=gSb+gSb+09ret00). The corresponding SFR is shown in Fig[TT](right panel). Each panel is 100 kpc x 100 kpc in size. 
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Fig. 14. Evolution with time of the gas lagrangian radii contain- 
ing 10% (black), 25% (red), 50% (green) and 75% (blue) of 
the gas mass, in units of the corresponding gas radii of the iso- 
lated galaxies. Left panel: local gSbgSb09ret00 interaction; right 
panel: gas-rich gSb+gSb+09ret00 encounter. In both panels, the 
lagrangian radii of each of the two galaxies involved in the inter- 
action are shown (solid and dashed lines, respectively). 



right panels in FigfTTIandfTSli. In turn, the amplitude of the first 
SFR peak, just after the first close passage between the two 
systems, is higher in the case of gas-rich interactions than for 
their gas-poor counterparts. These trends are common to all the 
simulations performed with the Tree-SPH code. 



3.5.1. Disk fragmentation 

In this section, we want to understand the mechanism respon- 
sible of stimulating higher levels of star formation at the first 
pericenter passage in gas-rich systems rather than in local ones. 
The first hypothesis that we have tested is that the high star for- 
mation rates found in gas-rich systems could be due to a more 
pronounced gas inflow in the central regions. In order to under- 
stand if this is the mechanism at work, we evaluated the radii 
r\o, r25, r^o and rvj containing, respectively, 10%, 25%, 50% and 
75% of the gas mass present in the local interacting system 
with id=gSbgSb09ret00 and in the coiTesponding gas-rich pairs 
gSb+gSb-i-09ret00. The evolution with time of these radii is 
shown in Fig[T4l As usual, in order to separate the effects due 
to the tidal interaction from secular evolution, the radii are nor- 
malized to that of the corresponding isolated systems. In both 
cases (local and gas-rich systems), it is evident that before the 
first pericenter passage, all these values are unchanged with re- 
spect to the isolated systems (indicating that at these times, the 
interaction is still not efficient in perturbing the gas dynamics). 
But then, just after the pericenter passage (t=400 Myr), the disks 
of the two gSb galaxies and of their gas-rich gSb+ counterparts 
are quite perturbed: this gives rise to a compression of the inner 
radii (rio, r25) of a factor of ten in both cases (see, respectively, 
left and right panel in FigfT4l). But the gas is compressed by the 
same amount in both pairs. This means that the inflow could not 
be responsible, alone, of driving different star formation rates, 
relative to the isolated systems, at this time: indeed the same 
amount of compression for local and gas-rich systems should 
stimulate the same amount of star formation enhancement. 

The important difference between the local interaction and 
that involving gas-rich systems resides in the fact that these lat- 
ter, being more unstable, are more prone to fragment under the 
tidal effects of the companion galaxy. In this sense, for these 
simulations, we should talk of fragmentation-driven starbursts 
rather than inflow-driven ones. In Fig[T3] for example, it is clear 
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that the first close passage of the two galaxies is accompanied by 
the formation of many gas clumps in the disks, which are also 
site of intense star formation. This fragmentation is less evident 
for the gSb+ galaxies evolved in isolation, thus indicating that 
this disk instability is largely a result of the encounter. Note also 
that the formation of gas clumps in the galaxy disks was less ev- 
ident for local interacting pairs (see maps in DM07). As we dis- 
cuss in more details in Appendix ICl in the case of some isolated 
gas-rich galaxies, the more important is the disk fragmentation, 
the higher is the star formation rate in the galaxy. This indicates 
that the disk fragmentation can strongly contribute to the great 
SFR enchancement found for gas-rich spirals at the pericenter 
passage. 

3.5.2. Effects of feedback 

If, on the one side, gas-rich systems present higher star formation 
enhancements than their local counterparts at the pericenter pas- 
sage, on the other side, in the merging phase, gas-rich galaxies 
show in general lower star formation rate^ than local interacting 
pairs (see FigsfTTIandfTZb. 

What is the physical mechanism that can inhibit a starburst in the 
merging phase of these gas-rich gSb+ spirals? As we checked, 
before the coalescence of the two systems, the gas amount is still 
enough to produce a considerable enhancement of the star for- 
mation rate. So, why is this gas reservoir not sufficient to drive 
important bursts of star formation? To answer this question, first 
of all we analyzed the evolution of the gas component during 
the gSb+gSb+09ret00 merger FigfT3] shows some maps of gas 
and star forming regions during the interaction. Since the spin of 
the two galaxies is antiparallel to the orbital angular momentum 
(retrograde interaction), the morphology of the two disks is only 
sligthly affected after their first close passage: no long tidal tails 
are formed, only a bridge connects the two systems. Note also 
that at that epoch (t=600 Myr), just after the first close passage, 
the SFR reaches its maximum value. As it can be seen from the 
maps, in this phase, the star formation is quite extended along 
the two disks, from the central regions to spiral arms. As the two 
galaxies begin to approach one another once again, the gaseous 
component tends to collapse in the central disk regions (t=1.3 
Gyr - 1 .4 Gyr), thus giving rise to the SFR enhancement found 
in FigHT] But just after this strong inflow, the gas expands, re- 
ducing drastically the star formation rate. 

The response of the gaseous component after the strong in- 
flow is a consequence of the kinetic feedback released by super- 
novae explosions. As described in DM 07, energy injection into 
the I SM is included following the work by M ihos & Hernquisti 
( 11994b.) . by applying a radial kick to velocities of neighboring 
gas particles. In DM07, performing a set of simulations of iso- 
lated galaxies, we found that a rapid growth of the vertical thick- 
ness of the gaseous disk can be prevented if the total amount of 
kinetic energy received by a gas particle, due to the contribu- 
tion from all neighbors, corresponds to less than 1 km/s. FigfTOl 
clearly shows that this is the case also for gas rich systems: in- 
deed, with such an amount of feedback, the vertical thickness of 
the gaseous disk remains quite constant during the evolution. 
Evidently, the amount of feedback acquired by a gas particle 
depends on the local star formation rate: particles that lie in a 
actively star forming region receive an energy injection greater 
than that acquired by particles in the disk outskirts. 
Galaxy interactions, perturbing the star formation evolution, also 
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Fig. 15. Probability to have a kinetic feedback larger than the 
value specified on the x-axis for gas particles, during the en- 
counter with id=gSb+gSb+09ret00 . Different curves corre- 
spond to different epochs during the interaction, as explained in 
the figure. 

change the amount of kinetic feedback received by the ISM. 
FigHSshows the probability of finding gas particles receiving an 
amount of feedback greater than a certain value, in an interval of 
values ranging from 10"^ to 6kms-i, for the gSb+gSb+09ret00 
encounter. In this plot, for each particle, the kinetic feedback has 
been evaluated as the total velocity impulse Av received by a gas 
particle, due to the contribution from all its neighbors. 
In the early phases of the interaction (t=100 Myr), when the two 
systems are sufficiently far away from each other to be consid- 
ered isolated, 50% of the gas particles receive a total radial ve- 
locity kick from neighbors which is less than 0.03 kms ', only 
1/20 having Av greater than 0.1 kms '. As the interaction pro- 
ceeds, and the two galaxies approch the pericenter (t=400 Myr), 
the amount of kinetic energy received increases. This trend con- 
tinues until t=600 Myr, when the peak of SFR occurs. But the 
strongest increase in the kinetic feedback of the system is ob- 
tained in the merging phase: at t=1.4 Gyr, indeed, 50% of the 
gas particles receive a velocity kick greater than 0. 1 kms"' . Note 
that such a high feedback for such a high fraction of gas particles 
is due to the strong inhomogeneity of the star forming regions. 
In other words, the strong gas inflow produces a highly concen- 
trated star forming region, where gas particles acquire conspic- 
uous radial kicks from neighbors. This rise in the kinetic energy 
determines the gas expansion found in Fig[T3]and the subsequent 
decrease in the star formation rate of the merger. 

4. Comparison with PIVI-SP simulations 

Are the main results found for Tree-SPH simulations still valid 
for PM-SP ones? Or, in other words, to what extent do the re- 
sults presented in the previous section depend on the numerical 
techniques adopted? Let us first recall the various assumptions 
that differ between the Tree-SPH and PM-SP models, and could 
a priori cause major differences in the star formation history: 

- A^-body and gas dynamics schemes are completely different 
(tree- versus grid-based). 

- The initial conditions differ: axisymmetrical in the Tree-SPH 
sample, with evolved spirals in the PM-SP models. Dark 
matter extents are also different. 

- A local Schmidt law was used in both models, but applied in 
a different way: in the Tree-SPH simulations, the gas density 
is computed using the adaptative tree structure, implicitely 
assuming a scale-free Schmidt law, while the gas density is 



14 



P. Di Matteo et al.: On the frequency, intensity and duration of starburst episodes triggered by galaxy interactions and mergers 



computed on a grid with constant resolution in the PM-SP 
simulations. 

- Different star formation schemes were used only in the PM- 
SP models. 

- One of the galaxies in the pair is always gas-free in the PM- 
SP models, while different morphologies (and gas contents) 
have been adopted for Tree-SPH simulations. 

Some other differences are not expected to influence much 
the star formation: 

- The orbital inclination is not varied in the PM-SP models, 
but according to the Tree-SPH models this parameter has lit- 
tle influence on the statistical properties of star formation 
(Figs.|6]to[8). 

- Star formation feedback is included only in the Tree-SPH 
models. 

We now discuss the robustness of the main results obtained 
from the Tree-SPH models by comparing them with the results 
of PM-SP simulations. Keeping in mind that these two sam- 
ples contain interactions with different orbital parameters (im- 
pact parameter, relative velocities, disk inclination), to compare 
them properly we have considered that encounters have differ- 
ent likelihoods, with a probability p cc b^V, b being the im- 
pact parameter of the encounter and V the relative velocity of 
the two systems. This choice is motivated since the collision rate 
is proportional to the relative velocity V and to the cross-section 
nb^. We have also tested other possibilities, including a depen- 
dence on the inclination / of the galaxy disks (for example, by 
means of a factor oc sin{i)), but this did not change significa - 
tively the results. We refer t he re ader torBournaud et alj (l2007bl) . 
iMartig & BournaudI (|2007|) and iMihosI (|2004') for a further dis- 
cussion on statistical weighting of interacting galaxy models. 
Results obtained this way, shown on Figs. [16] and [17] allow a 
statistical comparison of the two numerical datasets, as well as a 
statistical comparison of simulations with observations. 



4. 1 . Frequency and intensity of merger-induced starbursts 

Comparing the statistical distribution of the maximum star for- 
mation enhancement of local interacting pairs, it is striking to 
see how the results are in agreement, even adopting different 
numerical approaches. Indeed, Tree-SPH simulations and PM- 
SP ones with identical star formation recipes (Schmidt law, only 
density-dependent) have very similar distributions (see FigfTSb: 
in both cases, it results that interactions and mergers, in gen- 
eral, produce moderate enhancements in the star formation rate 
of the pairs, while strong starbursts are rare. Indeed, less than 
15% in each set of simulations show a maximum star formation 
rate greater than 5 times that of isolated systems. The median 
value of the SFR enhancement is 2.8 in the Tree-SPH dataset, 
and 3. 1 in the PM-SP dataset. 

Moreover, these results are not dependent on the models 
adopted to compute local star formation rates, as shown by 
the comparison of the various subset of PM-SP simulations. 
The l-sFR <^ ^gas^ model (compatible with observations, see 
Sect. |2.3| and Fig.[T]) indeed shows a statistical distribution quite 
similar to the Schmidt law model. Higher probabilities of strong 
starbursts are found with the Y^sfr ^ "^gas^'^^ model: as dis- 
cussed in Sect |2.3l (and Fig. [1]), this parametrization is barely 
compatible with observations and should be considered as an 
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Fig. 16. Probability to have a maximum SFR larger than the 
value specified on the x-axis for major mergers. The SFR is nor- 
malized to that of the corresponding isolated galaxies. The grey 
curves correspond to Tree-SPH simulations, the red curves to 
Particle mesh-sticky particles ones. Different lines correspond 
to simulations employing different recipes for modeling star for- 
mation: Y-sFR °^ ^eas^'^ (solid lines), I-sfr °^ '^gas^ (dashed line), 
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Fig. 17. Probability to have a star formation rate larger than two 
times the star formation of the isolated galaxies for a time larger 
than the value specified on the x-axis during major galaxy merg- 
ers. The grey curves correspond to Tree-SPH simulations, the 
red curves to Particle mesh-sticky particles ones. Different lines 
coiTespond to simulations employing different recipes for mod- 
eling star formation: Hsfr Sgo.,'''' (solid lines), 1,3 fr ^ga.?^ 
(dashed line), I-sfr ^gas^'^^ (dot-dashed line). 



upper limit to the actual starburst efficiency, hence confirming 
our earlier conclusions. Even in this favorable case, the major- 
ity of the encounters (60%) leads to a maximal star formation 
enhancement smaller than a factor of five. 

Starbursts stronger than a factor 5 compared to the isolated 
galaxies are rare for the most realistic star formation models 
(Schmidt law and "Lsfr ^gas^ model) in both the Tree-SPH 
and PM-SP simulations. However, in more detail, one can note 
a difference for very high efficiencies: some starbursts reach 
factors 10-20 in the Tree-SPH dataset, while no similar case 
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Fig. 18. Same as Figs. [16] and [17] for very gas-rich galaxies (gas 
fraction in the disk: 50%). The grey curves correspond to Tree- 
SPH simulations, the red curves to Particle mesh-sticky particles 
ones. Large star formation rates in Tree-SPH simulations are at- 
tributable to instabilities resulting in the collapse of massive star- 
forming clumps during the interaction (see text). The effect of 
this internal fragmentation was removed in PM-SP simulations 
so that star formation is driven only by the merger-induced grav- 
ity torques; starbursts are then much more modest. 

is found in the PM-SP dataset. An explanation might be that 
the PM-SP sample is smaller and has randomly missed these 
most efficient cases, however the sample should have been large 
enough to obtain a few cases in this range. This difference must 
then have a physical cause, which could be one of the following: 

- Some peculiar and very favorable conditions are not included 
in the PM-SP sample. In particular, perfectly coplanar en- 
counters between two gas-rich galaxies can enhance the star 
formation rate by direct collision between the two gas com- 
ponents when the two disks overlap: this will create high 
density regions and dissipate the gas angular momentum 
very rapidly, triggering the gas infall towards the center Such 
situations should, however, be rare and correspond only to 
very specific situations. 

- This could also result from the different application of the 
Schmidt law: in the Tree-SPH code, the Schmidt law is as- 
sumed to be scale-free and the gas density is computed on 
the tree structure. In the PM-SP code, the gas density is com- 
puted at a fixed scale on the Cartesian grid. Very dense small- 
scale structures will be resolved only in the first case, while 
their density will be smoothed at the grid scale in the second 
case. Which model is the most realistic depends on the phys- 
ical origin of the Schmidt law. If this la w comes from the 
turbulent density distribution in the ISM (' Elmegreenll2002l) 
it should be scale-free just like in the Tree-SPH code. If the 
Schmidt law r esults from the d ynamical timescale of gas 
clouds 1 / yjGp (iElmegreenl2002l) . it should not be scale free, 
and the density should be computed when clouds begin to 
form on rather large scale, not at smaller scales when they 
have begun to collapse. In this case, the Tree-SPH may over- 
estimate the star formation efficiency by resolving density 
increases at small scales. 

Some differences concerning the most efficient starbursts can 
thus be caused either by different numerical techniques, or by 
different physical assumptions on star formation. They never- 
theless affect only rare cases or have limited amplitudes. We 
note in particular that, for the moderate starbursts up to a factor 
^ 4 in the relative SFR (that are the most frequent), the PM- 
SP model finds higher average efficiencies than the Tree-SPH 
one (FigfTSIl. This could be caused by the absence of feedback 
in the PM-SP model, as is indeed expected that feedback should 
regulate the star formation activity. Moreover, the differences ex- 
pected by including feedback or not in the two models might be 



somewhat compensated for by some other assumptions: in the 
PM-SP model, not including feedback favours an increase of the 
star formation activity, but on the other hand evaluating the gas 
density on the fixed resolution grid smoothes the density peaks 
and should somewhat regulate star formation compared to the 
Tree-SPH method. 

Overall, the Tree-SPH and grid-sticky particle simulations 
agree that strong starbursts are not frequent in major galaxy 
mergers, with a median starburst efficiency of about a factor 3 
compared to isolated galaxies. This conclusion is robust versus 
numerical methods. It is also robust versus morphological types 
of merging galaxies (from the Tree-SPH sample) and the physi- 
cal assumptions on the model used to compute local star forma- 
tion rates (from the three PM-SP subsamples). 

4.2. Duration of merger-induced starbursts 

The duration of the starbursts triggered by galaxy mergers is 
found to be generally smaller than 500 Myr in the Tree-SPH 
models. The PM-SP models have a somewhat different statisti- 
cal distribution of the starburst duration (see Fig. [TtI i that also 
depends on the adopted star formation model. However, both 
codes and all star formation models are in agreement on the main 
conclusions: no more than 20% of interaction/merger-induced 
starbursts have durations longer than 500 Myr, the Schmidt law 
in PM-SP models giving even no such case. The Schmidt law as 
well as the I-sfr ^gas^ have median durations of 200 Myr, 
and average durations of 200-300 Myr. This is obtained for 
disk stellar masses of several 10'" Mq, thus durations should 
be even shorter for smaller disks or higher redshift disks that 
have shorter dynamical timescales. Also note that, as discussed 
in Sect. 13.31 the burst duration diminishes as the star formation 
enhancements increase, so that the values found in the left panel 
of Fig [17] should be taken as upper limits if they are compared 
to observations. Let us note that observations are biased towards 
the strongest starbursts, because of brightness selection effects. 

4.3. Tlie case of gas-ricli gaiaxies 

As discussed above, the Tree-SPH and PM-SP simulations show 
a good agreement when modeling local galaxy interactions. The 
gas fraction, varied from 10 to 30% in the Tree-SPH simulations, 
does not seem to have a major impact on the relative intensity 
of merger-induced starbursts: disks that contain more gas have 
higher SFRs both before and during mergers, and the ratio of 
enhancement is more or less unchanged (see also DM07). 

For particularly gas-rich systems (gSb-t- models, assumed 
to represent gas-rich high-redshift galaxies), the situation is 
much more complex, as shown in Fig[T8] The PM-SP mod- 
els show roughly unchanged distribution of the starburst rela- 
tive efficiency and duration, while the Tree-SPH models show 
an increase of both the efficiency and duration of starbursts. 
As we saw in Sect. 13.51 in the case of Tree-SPH simulations, 
it is mainly in the first phases of the interactions that gas-rich 
pairs show a higher star formation rate than galaxies with lower 
gas fractions, while in the late merging phase star formation 
activity is lowered when the gas fraction is increased, due to 
stronger feedback, or earlier gas consumption effects. We have 
also seen that, for gas-rich spirals, the star formation enhance- 
ment at the pericenter passage is mainly due to disk instabili- 
ties and clumps formation, rather than to stronger gas inflows. 
In other words, for gas-rich systems, it is a "merger-driven frag- 
mentation' rather than "merger-driven inflows" that drives the 
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Fig. 20. Gas (top) and stars (bottom) maps for a gSbO-gSb+ merger (with an impact parameter of 62.5 kpc and an encounter velocity 
of 150 km s"') for the PM-SP model. Each box is 70 kpc x 70 kpc in size. 
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Fig. 19. Face-on and edge-on views (left : stars, right : gas) of 
the initial conditions for the gSb+ gas-rich galaxy in the PM- 
SP models (each box is 30 kpc x 30 kpc in size). Note that 
gas clumps here ( i) already exist in the initial conditions while 
they form during the merger in the Tree-SPH simulations and ( ii) 
these clumps are smoother, less compact and less massive than 
those formed in the Tree-SPH models, owing to the reduced dis- 
sipation rate and non-adaptative resolution. 



star formation in the Tree-SPH simulations with 50% gas. The 
star formation enhancement driven by disk fragmentation is also 
responsible of the longest burst duration found in the right panel 
of Fig[T7](a more detailed discussion on this point can be found 
in Appendix). 



The different behaviour of PM-SP simulations is due to the 
fact that in order to prevent an excessive disk fragmentation in 
the galaxy evolved isolated, the gas dissipation in these simu- 
lations had been decreased (see Section 2.1). This essentially 
prevents the formation of dense gas clumps, sites of intense star 
formation. Actually, the gas-rich PM-SP models are shown in 
Fig. [19] for the starting conditions and Fig.|20]for a merger se- 
quence: some clumps of gas are also present, but owing to the re- 
duced dissipation they are less massive and less dense that those 
seen in the Tree-SPH simulations. Moreover, they are present in 
the starting conditions and not formed during the merger. There 
is no additional fragmentation during the merger, and the relative 
SFRs measured in these simulations result from merger-driven 
inflows. And then we find in these models that the relative SFR 
enhancement by mergers, as well as the duration of these bursts, 
are rather unchanged compared to galaxies with lower gas frac- 
tion. 

The case of gas-rich galaxies is then more difficult to study 
and different assumptions will lead to largely different evolution. 
However, when one disentangles what drives the star formation 
activity in each model, one can likely conclude that ( i) merger- 
driven gas inflows do not trigger more efficient starbursts in gas- 
rich merging galaxies compared to local galaxies, but ( ii) other 
processes like internal clump formation, which itself can be pro- 
voked by a galaxy interaction, can trigger stronger and longer ac- 
tivities of star formation. We nevertheless note from Appendix A 
that such gas clumps can also occur without mergers/interactions 
being needed, then triggering comparably high starbursts in non- 
interacting systems. Internal instabilities may then trigger star 
formation in all kinds of galaxies and not necessarily enhance 
the star formation activity in the interacting ones. 



5. Discussion 

The role played by interactions and mergers in affecting the 
evolution of galaxies has been widely studied in the last decades 
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by means of numerical simulations. The first numerical models 
of interacting galaxies in which all the dynamical components 
(bulge, disk and massive dark halo) were taken into account 
by means of fully N- b ody simulati ons ( with ~ 10'*) were 
developed by iBarnesI (1 19881 1 19921) and iHernquist & Spergell 
(Il992h . In the same years, the first attempt to incorporate the 
dynamics of the interstellar medium and star formation recipes 
into m odels of interactions between disk galaxies (Mihos et aL 
Il992h showed that the global SFR in merging galaxies could 
be enhanced by an order of magnitude with respect to pre- 
interaction levels, for several hundred Myrs, with star formation 
mostly occurring in the galaxy central regions. 

These results w ere c onfir med and extend ed by 

iMihos & HernquistI ( Il994ah and iMihos & HernquistI (Il996h . 
who incorporated s tar formation algorithms in a Tree-SPH code 
dMihos & Hemquis t 1994b) in order to investigate the trigger- 
ing of starbursts in merging disk galaxies. Their main findings 
showed that the properties of the merger-driven starburst depend 
on the morphology of the interacting systems, being galaxies 
without bulges more prone to disk instabilities in the first 
phases of the interaction, while central bulges tend to stabilize 
galaxy disks and to postpone gas inflows and the subsequent 
starburst to the final merging phase. These works also showed 
that the merger-driven starburts are very short, with a duration 
of approximately 50 Myr, but not rare (the SFRs of interacting 
systems being typically two orders of magnitude greater than 
that of the two corresponding isolated galaxies). While the sim- 
ulations discussed in this paper and in DM07 present the same 
dichotomy in the response of early-type and late-type galaxies 
to interactions-driven instabilities (see, for example. Figs. 10 and 
13 in DM07), contrary to these works, our findings show that 
the most intense star formation enhancements (a factor greater 
than 10 with respect to that of isolated systems) are not only 
short, but also r are. Lower star formation enhancements with 
respect to that o f lMihos & Hernquis3 (11996 ) were also found by 
ISpringell (l2000l) : even adopting a more sophisticated feedback 
with respect to that employed by iMihos & Hernquisi d 19961) 
and by us, they find a star formation rate at maximum ten times 
greater than that of isolated system s, so an order of magnitud e 
lower than the ones proposed by IMihos & Hernquisi (Il996h . 
The evolution of the star formatio n rate in the M ice, studied by 
means of numerical simulations bv lBarnesI ( |2004|) . also suggests 
star formation rates at most ten times greater than pre-interaction 
levels. Interestingly, this work compares different star formation 
prescriptions, showing how the extensions of the star forming 
regions can be affected when taking into account also gas shocks 
in the modeling. The different response of the star formation 
evolution in interacting galaxies with or without bulges was 
also confirmed by Springel et al. (2005). In this work, a stability 
analysis of isolated disk galaxies was performed, showing the 
difficulty to prevent disk instabilites when a large isothermal 
gas fraction is adopted. This is in accordance with the evolution 
of our gas-rich systems, which, in fact, shows that strong 
disk instabilities can occur unless gas dissipation is artificially 
reduced (as is done for gas-rich galaxies in PM-SP simulations). 
They also pointed out the role played by accreting black holes 
in quenching star formation during the final galaxy coalescence. 

A simple stellar feedback model was implemented by 
ICox et ai] (I2OO61) . in their simulations of disk galaxy major merg- 
ers. They adopted a feedback model which stores SNe energy 
within individual gas particles and dissipates it on a time-scale 
defined by t//, and according to an equation-of-state in star form- 
ing regions set by a n parameter The adopted values of r ft, and 
n determine different absolute star formation rates both for iso- 



lated and merging systems, but, for the most plausible feedback 
prescriptions (i.e. feedback not too low), the choice of different 
parameters does not strongly affect the relative star formation 
rates (see Table 3 in that paper). They also pointed out that the 
adoption of an isothermal equation of state for the gas can over- 
estimate the star formation rates of merging pairs, because of the 
underestimation of shock heatings in low gas density regions. 

In our simulations we did not focus on the dynamics of gas 
falling-back from tidal tails, or similarly the gas potentially ex- 
pelled by strong winds and falling-back at later-stages. Such pro- 
cesses could contribute to increase the star formation activity af- 
ter the merger, but the dynamical timescale in the outer regions 
is large so this will likely occur only after the merger remnant 
is well relaxed rather than during the merger itself. Also, these 
mechanisms would be expected to feed a large-scale disk rather 
than a dense central concentration, so the associated star forma- 
tion would last long with a moderate efficiency instead of con- 
sisting of a starburst activity. Dedicated studies would be needed 
to fully constrain the role of the late gas fall-back on merger 
remnants; in any case the associated star formation is not related 
to the merger-driven inflow process that we have studied in the 
present paper. 

6. Summary and comparison to observations 

6.1. At low redshift 

We have infeiTed from the two simulations dataset presented in 
this paper that, for galaxies with gas content typical of Local or 
low-redshift galaxies (typically 15% of gas in the disk and up 
to 30%), the starbursts induced by galaxy major mergers have a 
moderate intensity, star formation rates being rarely enhanced by 
factors larger than 5 compared to isolated galaxies, even at the 
peak of the starbursts. This result appears to be robust because it 
is found in both numerical methods used (Tree-SPH and PM-SP) 
and is not largely affected when different assumptions are made 
for the determination of the local star formation rate. 

A large-scale tidal field from dense cosmological 
structures could further enhance merger-induced starbursts 
(iMartig & Bournaudl2007l) . especially at high redshift. This was 
not accounted for here, but the average effect does not exceed 
a factor of two. It obviously does not affect at all mergers in 
the field or low-mass groups. Also, accounting for feedback 
or not leaves the statistical conclusions about unchanged when 
our two codes are compared. The Tree-SPH dataset shows that 
the star formation activity is comparable in mergers and fly-by 
interactions. We here focused mainly on int eractions of galaxies 
of comparable masses, but Cox et al.l (|M)8) show that minor 
interactions are much less efficient than major ones to trigger 
star formation. 

Strong starbursts can still result from galaxy major interac- 
tions and mergers. About fifteen percent on the galaxy interac- 
tions in our models have a merger-triggered starburst with a rela- 
tive efficiency higher than 5. To our knowledge, there is no other 
mechanism susceptible to trigger such intense starbursts in low- 
redshift disk galaxies, and this agrees well with the vast major- 
ity of LIRGs and ULIRGs in the local univer se sh owing signs 
of interactions or mergers dSanders & Mirabellll996t iDuc et aH 
1997). More generally, nearly all mergers in our model lead to 
a significant enhancement of star formation of at least 50%, and 
in the majority of the models the star formation rate is at least 
doubled during the merger: this also agrees well with the fact 
that obs erved mergers show a general trend to enhanced sta r for- 
mation dOeorgakakis et al.ll200d:ICharmandaris et al.ll200lh . 
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However, merger-induced starbursts have on average a lim- 
ited efficiency. Strong starbursts are rare among simulated merg- 
ers and interactions, and this again is in agreement with obser- 
vations: observed starbursts at low redshifts are merger-driven, 
but mergi ng galaxies do not necessarily show high star forma- 
tion rates dBergvall et al.ll2004l:lLi et al.|[2007l) . Independently of 
the code used and star formation recipe adopted, w e find a me- 
dian v alue for the maximal relative SFR of about 3. I jogee et al.l 
(I2007h have studied star formation in interacting galaxies in 
GEMS at low and moderate redshifts (z ~ 0.24 to 0.80) showing 
that the average SFR of strongly distorted interacting/merging 
massive galaxies is only modestly enhanced with respect to nor- 
mal undisturbed systems, in agreement with the results from our 
models. Note that simulations can result in a slightly higher SFR 
enhancement than observations because the instant of the maxi- 
mal SFR is selected in simulated pairs. 

6.2. At high redshift 

At higher redshift, in particularly gas-rich galaxies, the situation 
is more complex. It appears that merger-induced gas inflows are 
not more efficient trigger of star formation than they are at low 
redshift. The star formation activity can be enhanced in much 
larger proportions, and over much larger durations, if the gas- 
rich disks become Jeans unstable and form dense gas clumps 
during the interaction, as is the case in the Tree-SPH simula- 
tions. However, comparison with the PM-SP simulations and 
other simulations presented in the Appendix indicate that this 
is a rather particular situation. Such clumps do not form in disks 
that are too stable - the starbursts have low average efficiencies 
- while in disks that are too unstable the star-forming clumps 
form even without the interaction - then only a moderate addi- 
tional starburst can be triggered if an interaction occurs. 

Under specific condition£3 where disks are at the stabil- 
ity limit before merger occurs, merger-driven fragmentation can 
lead to strong merger-driven starbursts. However, there are evi- 
dences that disks can become clumpy as a result of their inter- 
nal evolution only, without mergers/interactions being required 
for that (Elmeg reen et al. 2007; Bournaud et al. 2007a) (see also 
iBournaud et al] ((2008)), which can then be another way to trig- 
ger star formation without galaxy interactions. 

Observationally, the role of mergers in the star formation 
history at high redshift is still debated (see Introduction). There 
are however severa l studies (e.g. Bell et al. 2005) who find that 
this role is limited. iDaddi et"aiT(l2008l) find long duty cycles for 
ULIRGs at high redshift, which are longer than the typical dura- 
tion of merger-induced starbursts in our simulations, suggesting 
that another mechanism rather than interactions/mergers is the 
main trigger of starbursts and AGNs there. Large star formation 
rates could then be triggered by internal evolution, or simply re- 
sult from the presence of large gas reservoirs around the most 
actively star-forming galaxies CDaddi et al.,2008.) . 

7. Conclusion 

In this work we have analyzed the relation between galaxy 
interactions and star formation enhancement by means of a vast 
number (~ 1000, in total) of simulations, varying the numerical 
code adopted as well as the numerical recipes to model star 

Note, however, that cold gas accretion from filaments in the near 
environment of galaxies has not been taken into account in this work. 
We plan to study this aspect in future works, as it can affect the galaxy 
disk stability, as well as its star formation rate. 



formation. The main results of this study are the following. 

- At low redshifts, interactions and mergers, in general, pro- 
duce moderate enhancements in the star formation rate of 
the pairs, while strong starbursts are rare. This result does 
not depend either on the numerical code adopted, or on the 
recipes used to account for star formation. For star formation 
prescriptions compatible with observations (1,sfr ^gas^'^ 
or HsFR "^gas^), the majority of the encounters (~ 85%) 
leads to a maximal star formation enhancement less than a 
factor five. As discussed in Sect l6.1l these res ults are in good 
agreement with several observational s tudies dBergvall et al.l 
l200l:llTetani2007tlJogee et al.ll2007h . 

- The duration of the moderate starbursts is generally smaller 
than 500 Myr for Tree-SPH models. Even if PM-SP simula- 
tions have a somewhat different statistical distribution, both 
models show that no more than 15 % of interaction-induced 
starbursts have a duration greater than 500 Myr This would 
suggest that another mechanism has a major role in trigger- 
ing the acti vity of high-redshif t ULIRGs, which has a longer 
duty cycle dPaddi et al.ll2007al) . 

- Moving to higher redhifts (i.e. essentially increasing the 
gas mass fraction), interacting systems have absolute SFRs 
greater than local interacting pairs. 

- Inflow-induced starbursts during interactions and mergers 
are neither strongej^ nor longer than their local counter- 
parts. In turn, Jeans instability in gas-rich di sks can cause the 
formation of massive clumps, as ob served dElmegreen et al.l 
120071: IBournaud et al.ll2007all2008l) . The disk fragmentation 
can appear either in interacting either in isolated systems 
and the induced fragmentation-driven starbursts are charac- 
terized by high star formation enhancements as well as long 
bursts duration, but then interactions and mergers are not 
necessarily required to trigger this kind of activity. 

Additional numerical work is required to firmly establish 
the nature of high-z ULIRGs and understand the role of inter- 
nal factors compared to interactions and mergers in the trigger- 
ing of the star formation activity and/or AGN. Nevertheless, this 
work presents strong results in terms of frequency and duration 
of starbursts episodes triggered by galaxy interactions. Strong 
starbursts can occur in some major mergers, and this is con- 
siste nt with submillimit er galaxies being mostly major mergers 
(e.g. iTacconi et al.ll2008i) . but, on average, the classical merger- 
driven gas inflows do not seem to be an efficient process, and 
other mechanisms could have an important role too, like the frag- 
mentation of gas-rich disks, or the presence of lar ge gas reser- 
voirs as recently suggested by iDaddi et aTl d2008l) . These con- 
clusions hold well with observations showing that a large num- 
ber of actively star- forming galaxies are massive disks rather 
than mergers (Genzel et al. 2006; Forster Schreiber et al. 20061 
, Daddi e t al. 2007a; Shapi ro et al. 2008). El megreen et al. (200% 
have shown that resolved star-forming galaxies in the Hubble 
Ultra-Deep field at redshift 1 and above are largely dominated by 
primordial disks. These high- redshift disks outnumber merger 
candidates; they form stars at high rates because they are mas- 
sive and gas-rich: pervasive clumpy structures attest of the high 
gas content of these massive star- forming disks. These results 
certainly also help to disentangle the nature of massive star for- 
mation at high redshift, in the light of the limited role of mergers 
in systematically inducing extreme star formation rates, as sug- 
gested by our models. 

" The SFR being normalized to that of isolated galaxies. 
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Appendix A: Tree-SPH and PM-SP codes: some 
validity tests 

A.1. Tree-SPH 

The main parts of the Tree-SPH code used to run most of 
the simulations presented in this paper has been described in 
Semelin & Combes (2002). In particular, we employed the same 
evaluation of the gravitational forces, the same implementation 
of the SPH technique and of the star formation modeling as the 
one described in the above ci ted paper. Some validity t ests of this 
code have been presented in 'Se melin & Combes! (120021) . in par- 
ticular that concerning the collapse of an initially static, isother- 
mal sphere of self-gra vitating gas -a standard validation test for 
this type of codes, see lEvrardld 19881) and also lHernquist & Kat j 
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( 119891) : iThacker et al] (l2000h : ISpringel et all ( l2001h '). Here we 
present some other tests, about the dependency of the star for- 
mation rates on the numerical parameters adopted (integration 
time-step, number of neighbors A^, used in evaluating SPH quan- 
tities, value of the gravitational smoothing length e, and number 
of total particles in the systems). In particular, we have chosen 
one of the coplanar merger simulations described in Sect |3.5l (the 
one with id=gSb-i-gSb+09ret00) and we have performed four ad- 
ditional runs, changing: 

- the value of the integration time-step, taken as one half of 
that of the reference case, i.e. Af - 25Myr, 

- the value of the number of neighbors A^, adopted to evalu- 
ate the SPH quantities, performing one run with Ns = 30, 
instead of A^, = 15 : 

- the value of the gravitational smoothing length, assuming 
e - 350/7C instead of e = 280pc : 

- the total number of particles of the simulation, performing a 
high-resolution run with a total number of particles equal to 
N,oi = 960000, i.e. four times higher than the reference case 
value. 
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The results of this study are presented in Fig lA.ll left panel, 
where the relative SFR for the reference case (black curve) and 
for the additional tests are presented. As it can be seen, the two 
peaks of the SFR (corresponding to the time of pericenter pas- 
sage and to the merging phase) are only slightly affected by the 
choice of different numerical parameters, being the differences 
in the relative SFR values of ~ 10% at the first passage, and of 
~ 20% at most in the final coalescence phase. 

A.2. PM-SP 

Some validity tests have been run also for the PM-SP simula- 
tions. In particular, we have analyzed the dependency of the rel- 
ative SFR on the integration time-step, choosing six representa- 
tive mergers, which have a maximum relative SFR equal to 4.66 
on average. Dividing the initially adopted time step Af - IMyr 
by two, we have checked that the evolution of the relative SFR is 
the same (cfr Fig lA.ll right panel), with a maximum SFR which 
can be at most ~ 20% greater or smaller than the reference value, 
but with no systematic effects with respect to it, which seems to 
be a sampling effect rather than a physical difference. 

Appendix B: Some evolutions of absolute star 
formation rates during interactions and mergers 

In order to separate the contribution of tidal effects from secular 
evolution in determining the star formation history of the 
interacting pairs, all the star formation rates presented in the 
paper have been shown relative to that of the two corresponding 
isolated galaxies. Here we want to show the evolution of the 
absolute star formation rates during some encounters, in order 
to give some indications about the absolute rates involved. 
In Figs. IB.ll and IB. 2| this is done for some coplanar encounters 
involving two gas-rich spirals and two local Sb spirals, respec- 
tively (cfr Figs. [TT] and [12] for the evolution of their relative 
SFRs). In the case of the gSb-i-gSb-nOlretOO interaction, the 
merging phase occurs soon after the first pericenter passage 
and the burst of star formation {SFR ~ OOMo/yr) takes place 
between this two phases. At the same epoch (f ~ 450^^3^), the 
isolated gSb+ galaxy sustains a SFR of the order of 4 M^lyr. 
This leads to a normalized star formation rate (relative to that of 
the two gSb-H galaxies involved in the interaction) of about 15 
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Fig. A.l. Left panel: validity tests for the Tree-SPH simulations. 
The relative SFR evolution during an interaction is shown, either 
when using the standard values described in Sect l2.2.ll (solid 
thick black curve), or an integration time-step half of that of the 
reference case (solid thin red curve), or a gravitational smooth- 
ing length 1 .4 times higher (dashed green curve), or a number of 
neighbors A^, two times higher than that of the reference case 
(dot-dashed blue curve), or a number of total particles in the 
system four times higher than the reference value (dotted vio- 
let curve). Right panel: validity test for the PM-SP simulations. 
The relative SFR evolution during an interaction is shown, when 
using either the standard values described in 12.2.21 (solid thick 
black curve) or an integration time-step half of that of the refer- 
ence case (solid thin red curve). 



(see left panel in Figs IB. It . After this starburst, the SFR of the 
merger decreases rapidly, to the consumption of most of the gas 
available. 

The evolution of the star formation history of the 
gSb+gSb+09ret00 interaction is shown in the right panel 
of Fig lB.ll in this case the star formation rate of the interacting 
galaxies begins to increase rapidly after the first pericenter 
passage, peaking at about 36 M^lyr 200 Myr after the first 
encounter; a second burst of the same absolute magnitude 
occurs when the coalescence of the two galaxies takes place. 
As discussed in Sect l3.5l the amplitude of this second burst of 
star formation is limited by the strong kinetic feedback due to 
SNe explosions, which determines a rapid expansion of the gas 
component. 

A similar behaviour to that of gSb+ interactions is found for the 
SFR evolution of local galaxies (Fig. |B.2| l, except for the fact 
that these systems contain a smaller gas fraction (cfr Table ^ 
and so, ultimately, they show lower star formation rates during 
the interaction. 

Finally, the evolution of the absolute star formation rate of a 
representative interacting gSb-Sb pair and of the corresponding 
isolated galaxy in the case of PM-SP simulations is given in 
Fig[ 



Appendix C: Disk fragmentation and star formation 
enhancement 



In Sect |3.5l we found that gas-rich spirals in interaction are 
more prone to disk fragmentation than their local counterparts 
and that the induced formation of gas clumps is mainly respon- 
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Fig.B.l. Absolute star formation rate, versus time, for some 
coplanar mergers involving two gas-rich gSb+ galaxies (solid 
lines). The absolute star formation rate of the isolated gSb-H 
galaxy is also shown (dotted line). In each panel, the black ar- 
rows indicate, respectively, the first pericenter passage between 
the two galaxies and the merger epoch. See Fig[TT]for the evolu- 
tion with time of the relative SFR for these two encounters. 
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Fig. B.2. Absolute star formation rate, versus time, for some 
coplanar mergers involving two local gSb galaxies (solid lines). 
The absolute star formation rate of the isolated gSb galaxy is 
also shown (dotted line). In each panel, the black arrows indi- 
cate, respectively, the first pericenter passage between the two 
galaxies and the merger epoch. See Fig[T2]for the evolution with 
time of the relative SFR for these two encounters. 



sible for the star formation enhancement found with respect 
to isolated systems. To better clarify the role played by disk 
fragmentation in stimulating star formation, here we want to 
compare the star formation evolution of the gas-rich gSb-H 
galaxy (see FiglTO]for the gas maps) with that of two gas-rich 
systems (hereafter called gSb+ul and gSb+u2), having the same 
morphological parameters, but initial Toomre Q parameters 
lower than the gSb-n one. In more detail, the Toomre parameter 
Q, equal to unity for the gas component of the gSb-H galaxy, has 
been lowered io Q - 0.3 for the gSb+ul galaxy and to 2 = 0.1 
for the gSbH-u2 one. 

The effect of reducing the local disk stability on the galaxy 
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Fig. B.3. Absolute star formation rate, versus time, for a PM-SP 
merger (solid line). The absolute star formation rate of the iso- 
lated gSb galaxy is also shown (dashed line). The black arrows 
indicate, respectively, the first pericenter passage between the 
two galaxies and the merger epoch. 

evolution is shown in Fig lC.ll and IC.2I with respect to the 
stable case (FigfTOt. the gas disk fragments in many clumps, 
particularly evident for the gSb-i-u2 system. These clumps form, 
migrate towards the galaxy centers and, in som e cases, dissolve 
under the influence of the tidal field of the disk jBournaud et alj 
I2007ah . 

It is natural to ask if a higher level of gas clumps formation 
in the disk also results in a different star formation evolution for 
these isolated systems. As shown in Fig |C.3l this is indeed the 
case. Comparing the star formation rate of the unstable gSb+ul 
galaxy with that of the more stable gSb-n one, we found that dur- 
ing 3 Gyr of evolution, the unstable disk shows a SFR higher 
than that of the stable counterpart. This is especially true after 
t=2 Gyr, when some clumps in the gSb-i-ul galaxy are migrating 
in the central region of the system: in this phase, the ratio be- 
tween the two star formation rates increases up to a factor of 2.5. 
Even more striking is the gSbH-u2 case, the most unstable one: 
in this case, a lot of clumps form in the first Gyr of evolution, 
giving rise to a star formation burst which peaks at 6 times that 
of the stable gSb+ galaxy. 

The presence of gas clumps in the disk also affects the duration 
of the star formation enhancement, in the sense that the most 
clumpy is the disk, the longest is the duration of the star forma- 
tion enhancement. As shown in Fig |C.4| in general, the unstable 
galaxies sustain star formation enhancements greater than a cer- 
tain threshold (the threshold being varied from 2 M^/yr to 21 
Ma/yr) for a duration longer than the one of the corresponding 
stable galaxy. These few examples seem to show that the most 
the galaxy fragments under the effects of local disk instabilities, 
the greatest is the burst enhancement and the longest is the burst 
duration. 

Appendix D: Statistics for Tree-SPI-i simulations 

In this Appendix, some more details about the statistical distri- 
butions of the intensity, duration and integrated star formation 
rate of interacting galaxy pairs are given. In particular, three ta- 
bles are shown. The first one (Table ID. Il l presents the statisti- 
cal distribution of the maximum star formation rates for merger 
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Fig. C.l. Gas maps for the unstable gSb+ul galaxy, evolved in 
isolation. From left to right, maps are shown from t=500 Myr to 
t=3 Gyr, every 500 Myr. Both xy projection (bottom panels) and 
xz projection (top panels) are shown. Each box is 40 kpc x 40 
kpc in size. 



Fig. C.l. Gas maps for the unstable gSb+u2 galaxy, evolved in 
isolation. From left to right, maps are shown from t=500 Myr to 
t=3 Gyr, every 500 Myr. Both xy projection (bottom panels) and 
xz projection (top panels) are shown. Each box is 40 kpc x 40 
kpc in size. 
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Fig. C.3. Evolution with time of the star formation rates of the 
gSb+ul (dotted line) and gSb+u2 (dashed line) galaxies, nor- 
malized to that of the stable gSb-H system. 

and flybys, whose corresponding histograms are shown in Fig|6] 
The duration of the enhanced star formation rate is given in 
Table ID. 2l for two different SFR thresholds. We refer the reader 
to Sect l3.3l and Fig|7]for a discussion on this. Finally, in Table 
ID. 31 the quartiles of the distribution of the integrated star forma- 
tion rate, for mergers and flybys, are given. The corresponding 
histograms are shown in Fig (8] 
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Fig. C.4. Duration of a star formation enhancement greater than 
the value specified on the x-axis for the stable gSb-H galaxy (solid 
curve) and for the unstable gSb-nul (dotted curve) and gSb-i-u2 
(dashed curve) systems. 

Table D.l. Quartiles of the distribution of the maximum star for- 
mation rate (relative to the isolated case) for mergers and flybys. 
The first (^"25), second (P50) and third (^"75) quartiles are given, 
for all the different disk inclinations. See Fig|6]for the histograms 
of the distributions. 
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Table D.2. Quartiles of the distribution of the duration (in Myr) 
of enhanced SFR for the whole sample of interacting galaxies 
(mergers and flybys). Two thresholds are shown: relative SFR > 
2 and relative SFR > 5. The first (P25), second {P50) and third 
(f 75) quartiles are given, for all the different disk inclinations. 
See Fig|2]for the histograms of the distributions. 
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Table D.3. Quartiles of the distribution of the integrated star for- 
mation rate (relative to the isolated case) for mergers and flybys. 
The first {P25), second (Pso) and third (^'75) quartiles are given, 
for all the different disk inclinations. See Figl8]for the histograms 
of the distributions. 
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